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We study a strong suppression of the relative production rate {d — Au) / ijp—p) for inclusive high-pr 
hadrons of different species at large forward rapidities (large Feynman xp)- The model predictions 
calculated in the light-cone dipole approach are in a good agreement with the recent measurements 
by the BRAHMS and STAR Collaborations at the BNL Relativistic Heavy Ion Collider. We predict 
a similar suppression at large pt and large xp also at lower energies, where no effect of coherence is 
possible. It allows to exclude the saturation models or the models based on Color Glass Condensate 
from interpretation of nuclear effects. 

PACS numbers: 24.85.+p, 12.40.Gg, 25.40.Ve, 25.80.Ls 



I. INTRODUCTION 

In the proton(deuteron)-nucleus collisions, investi- 
gated at the Relativistic Heavy Ion Collider (RHIC), re- 
cent measurements of high-py particle spectra at large 
forward rapidities performed recently by the BRAHMS 
[1, 2] and STAR 3 Collaborations show a strong nuclear 
suppression. The basic explanation for such an effect 
has been based on an idea that in this kinematic region 
corresponding to the beam fragmentation region at large 
Feynman xp one can reach the smallest values of the 
light-front momentum fraction variable X2 in nuclei. It 
allows to access the strongest coherence effects such as 
those associated with shadowing or the Color Glass Con- 
densate (CGC). 

It was shown in [J, Q that a considerable nuclear sup- 
pression for any reaction at large xp (small X2) is caused 
by another effects, which can be easily misinterpreted 
as coherence. Such a suppression, for example, can be 
treated in terms of the Sudakov form factor reflecting the 
energy conservation. It is governed by the probability to 
produce no particles at large xp ^ ^■ 

Nuclear suppression at large xp can be also inter- 
preted alternatively, as a consequence of a reduced sur- 
vival probability for large rapidity gap (LRG) processes 
in nuclei, an enhanced resolution of higher Fock states 
by nuclei, or an effective energy loss that rises linearly 
with energy. It was demonstrated in refs. ^, ^ that it is 
a leading twist effect, violating QCD factorization. 

The BRAHMS Collaboration [ij in 2004 for the first 
time found a substantial nuclear suppression for high- 
Pt negative hadrons produced at large pseudorapidity 
ry = 3.2 (see Fig. [iJ. Because the data cover rather small 
X2 ~ 10~^, the interpretation of such a suppression has 
been tempted to be as a result of saturation [6, % or the 
CGC expected in some models |9j. 

Alternative interpretation of the nuclear effects occur- 



ring in the BRAHMS data [l] is based on energy con- 
servation implemented into multiple soft rescatterings of 
the projectile quark in nuclear matter [1, [Hj- Moreover, 
new data for neutral pions from the STAR Collaboration 
have been recently appeared at the same cm. energy, 
= 200 GeV but at still larger pseudorapidity 77 = 4.0 
demonstrating a huge nuclear suppression, which is more 
than a factor of 2 larger than at 77 = 3.2 manifested by the 
BRAHMS data. Although a minor part of this difference 
can be explained by the isospin effects, i.e. by a differ- 
ence in production of h~ and tt*' particles on deuteron 
target, there is still a large room for investigation of this 
huge suppression. It provides a good possibility to test 
our interpretation of such an effect. This represents one 
of the main goal of the present paper. 

Another very interesting effect following from our in- 
terpretation of nuclear effects at large x p and supported 
by available data is the xp scaling of nuclear suppression. 
It is in contradiction with X2 scaling, which is expected 
if the scaling represents the net effect of quantum co- 
herence. The detailed analysis of this xp scaling can be 
found in [1, Q collecting the data for production of differ- 
ent species of leading hadrons with small pt vap—A colli- 
sions at different energies covering the laboratory energy 
range from 70 GeV to 400 GeV. Similar x p scaling of the 
nuclear effects is observed from data on J/^E* and pro- 
duction measured by the E866 Collaboration at Fermilab 
[lOj compared to lower-energy data Besides, recent 
measurements of nuclear suppression for J/^P production 
in d- Au collisions by the PHENIX Collaboration [l2| 
at RHIC are consistent with xp scaling and exhibit a 
dramatic violation of X2 scaling when compared with the 
E866 data [13] ■ 

According to a;_F scaling [3, we expect similar nu- 
clear effects at forward rapidities also at lower energies 
when the onset of coherence effects is much weaker. It 
gives much less room for explanation of a strong nu- 
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clear suppression in terms of CGC. Because new data 
of high-pr hadron production at forward rapidities are 
expected also at smaller cm. energies ^ = 130 GeV 
and 62.4 GeV, corresponding predictions for nuclear ef- 
fects at large pT will play an important role for further 
verification of various phenomenological models based on 
the CGC. 

Besides, new data from the BRAHMS Collaboration 
for different species of high-pT hadrons produced at 
rj = 3.0 have been recently appeared. It allows to provide 
another probe of our model in investigation of significant 
nuclear suppression at forward rapidities. This coincides 
with the further goal of the present paper. 

The paper is organized as follows. In Sect. |TT] we 
present shortly a formulation of nuclear suppression in 
terms of Sudakov suppression factors adopted for mul- 
tiple parton interactions with the nucleus. Here we also 
mention about three different mechanisms of high-p^ par- 
ticle production. In the next Sect. IIIll we calculate the 
predictions for the ratio of particle production rate in 
d — A and p ~ p collisions as a function of px at large 
forward rapidities corresponding to BRAHMS and STAR 
experiments at RHIC. We analyze much stronger nuclear 
suppression obtained recently by the STAR Collabora- 
tion 3 at = 4.0 in comparison with the well known 
results from the BRAHMS experiment [l[ at ry = 3.2. 
We demonstrate that the model calculations without any 
free parameter are in a good agreement with data from 
the both collaborations. Finally we perform predictions 
for nuclear suppression at large forward rapidities also at 
lower energies, where no effect of coherence is possible. It 
allows to exclude the models based on CGC from expla- 
nation of a strong nuclear suppression. We also demon- 
strate an approximate exp(?7)/\/s {xp) scaling of nuclear 
effects in the energy and pseudorapidity range accessible 
by the BRAHMS and STAR Collaborations. The results 
of the paper are summarized and discussed in Sect. IIVI 



II. HIGH-pT HADRON PRODUCTION AT 
FORWARD RAPIDITIES: SUDAKOV 
SUPPRESSION, PRODUCTION CROSS 
SECTION 

Any hard reaction in the limit xp 1 can be treated 
as LRG process, where gluon radiation is forbidden by en- 
ergy conservation. If a large-x^? particle is produced, the 
rapidity interval to be kept empty is Ay = — ln(l — xp). 
Assuming as usual an uncorrelated Poisson distribution 
for radiated gluons, the Sudakov suppression factor, i.e. 
the probability to have a rapidity gap Ay, was developed 
in ref. 'i'l and has a very simple form. 



S{xp) = 1 — Xp 



(1) 



Nuclear suppression at xp 1 can be formulated as a 
survival probability of the LRG in multiple interactions 
with the nucleus. Every additional inelastic interaction 



contributes an extra suppression factor S{xp). The prob- 
ability of an n-fold inelastic collision is related to the 
Glauber model coefficients via the Abramovsky-Gribov- 
Kancheli ( AGK) cutting rules [l3| • Correspondingly, the 
survival probability at impact parameter b reads, 



hA 
LB.G 



exp[— ( 
5Z TT] 



hN 



TAib)] 



(2) 



where TA{b) is the nuclear thickness function. 

Assuming large values of hadron transverse momenta, 
the cross section of hadron production in dA{pp) colli- 
sions is given by a convolution of the distribution function 
for the projectile valence quark with the quark scattering 
cross section and the fragmentation function, 



d^a 
d^PT dr] 

d^qAip)] 



E 



d'^QT drj 



qT=PT/z 



dz fq/d(p){xi,qT) 



(3) 



where xi = qt e"^ / \/s- The quark distribution func- 
tions in the nucleon have the form using the lowest order 
parametrization of Cluck, Reya and Vogt [l^]. For frag- 
mentation functions we use parametrization from [l5l | . 

The main source of suppression at large px concerns 
to multiple quark rescatterings in nuclear matter. The 
quark distribution in the nucleus can be defined perform- 
ing summation over multiple interactions and using the 
probability of an n-fold inelastic collision related to the 
Glauber model coefficients with Gribov's corrections via 
AGK cutting rules ^1^. It has the following form: 



fq/Nix, 9t> M) = E ^) fq/Ni^^ It) 

n=0 

where the coefficients reads 

CTe// T{b, Z) 

Vnib,z) = 



(4) 



1 + (^eff T{b, Z) 



n+1 



(5) 



The effective cross section Ueff was evaluated in 

The valence quark distribution functions f^ij^ix, q^) in 
Eq. ([U are also given by the GRV parametrization [T^ 
but contain extra suppression factors, S{x)^ = (1 — x)^ 
([1]), corresponding to an n-fold inelastic collision, 



Qr,{x,ql) = Cnfq/N{x,ql)S{xy\ 



(6) 



where the normalization factors C„ are fixed by the Got- 
tfried sum rule. 

The cross section of quark scattering on the target 
da[qA{p)] / d^ qpdri in Eq. © is calculated in the light- 
cone dipole approach |16l. Il7|. Performing calculations. 
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we separate the contributions characterized by different 
initial transverse momenta of the projectile partons and 
sum over different mechanisms of high-pT production. 

Quark-diquark break up of the proton. Here 
we consider proton breakup remaining the diquark in- 
tact, p qq + q. We treat the diquark {qq} as point- 
like and integrate over its momentum. The correspond- 
ing kT distribution of the projectile valence quark, after 
propagation through the nucleus at impact parameter b, 
is is calculated using the dipole technique developed in 
refs. [H, [l^ (see also Q). This contribution dominates 
the low transverse momentum region fc^ < 1 GeV. 

Diquark break up qq — > qq- At larger kx the in- 
teraction resolves the diquark, so its break-up should be 
included. In this case the valence quark has much larger 
primordial transverse momentum. Its contribution is cal- 
culated in accordance with [l^, [H (see also Q). 

Hard gluon radiation q — > Gq. At large kx the 
dipole approach should recover the parton model [20j . 
where high momentum transfer processes occur (in lead- 
ing order) as binary collisions with the transverse mo- 
mentum of each final parton of order kx- Therefore, one 
should explicitly include in the dipole description radia- 
tion of a gluon with large transverse momentum that ap- 
proximately equilibrates fc^, i.e. the process qN — s- qGX. 
We employ the nonperturbative quark-gluon wave func- 
tion developed in |19l| .which corresponds to small gluonic 
spots in the nucleon [21, 22]. Details of calculation of the 
quark scattering cross section can be found in [4]. 



III. COMPARISON WITH DATA 

Several years ago the BRAHMS collaboration per- 
formed measurements of nuclear effects at pseudorapidity 
•q = 3.2 for production of negative hadrons with trans- 
verse momentum up to pT ~ 3.5 GeV. Instead of the 
usual Cronin enhancement a suppression was found, as 
one can see from Fig. [TJ 

Rather strong nuclear suppression of data at small pt 
has been analysed and interpreted in details in refs. [1, Q 
and such an analysis does not need to be repeated here. 
On the other hand we will concentrate on a study of 
nuclear effects at large pr- 

Note that the dominance of valence quarks in the pro- 
jectile proton leads to an isospin-biased ratio. Namely, 
high-pT negative hadrons close to the kinematic limit are 
produced mainly from d, rather than u, quarks. There- 
fore, more negative hadrons are produced by deuterons 
than by protons, and this causes an enhancement of the 
ratio plotted in Fig.[T]by a factor of 3/2. We take care of 
this by using proper fragmentation functions from p^ . 

Although the nuclear effects under discussion are not 
sensitive to pr dependence of the cross section for hadron 
production in p — p collisions, the model has been al- 
ready successfully confronted with p — p data from the 
BRAHMS experiment Q at r; = 3.2 in refs. 

As the next step, very important for verification of 
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FIG. 1: Ratio, Rd+Auipr), of negative particle production 
rates in d — Au and p — p collisions as function of pr at pseu- 
dorapidity r] — 3.2 vs. data from the BRAHMS Collaboration 
[^. Ratio Rd+Auipr) for neutral pion production at = 4.0 
vs. data from the STAR Collaboration Solid and dashed 
curves correspond to calculations with the diquark size 0.3 fm 
and 0.2 fm respectively. 



our model, we calculate nuclear effects employing the 
dipole formalism and the mechanisms described shortly 
in the previous section (see also Q). The results are com- 
pared with the BRAHMS data for the minimum-bias ra- 
tio Rd+AuiPr) [H in Fig- HI One can see that calculations 
are in a rather good agreement with data. 

In the same Fig. [T] we show also the STAR data for 
tt'^ production presented as the ratio Rd+AuiPr) at pseu- 
dorapidity ry = 4.0 3. Data demonstrate that suppres- 
sion is much stronger than at 7y = 3.2 observed by the 
BRAHMS Collaboration fl]. A part of this difference 
can be explained by an isospin-biased ratio followed from 
the dominance of valence quarks in the projectile pro- 
ton (deuteron). Whereas more negative hadrons are pro- 
duced by deuterons than by protons, for positive hadrons 
the situation is opposite. 

For production of negative hadrons, it leads to an en- 
hancement of the ratio Rd+AuiPr) by a factor of 3/2 in 
comparison with Rp^AuiPx)- However, for production of 
positive hadrons one arrives to a suppression by a factor 
of 3/4. As a result, for 7r° mesons one gets for Rd+AuiPr) 
a small overall suppression factor = 5/6, which is smaller 
than a factor of 3/2 for negative hadrons. However, such 
a difference following from the isospin effects can explain 
only a minor part from a huge difference in nuclear sup- 
pressions experimentally observed at 77 = 3.2 and 4.0 by 
the BRAHMS [H and STAR [| Collaborations, respec- 
tively. 

If one supposes to interpret the nuclear effects of high- 
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FIG. 2: Model predictions for nuclear attenuation factor 
Rd+Auipr) as a function of transverse momentum for pro- 
duction of tt" mesons at y/s = 200 GeV and different values 
of pseudorapidity "q from 3.0 to 4.0. For an illustration, full 
squares represent the data from the STAR Collaboration 0] 
obtained &t rj — 4.0. 



Pt hadron production at 77 = 3.2 in terms of CGC, such 
an interpretation should fail at larger 77 = 4.0, where 
observed suppression by the STAR Collaboration Q is 
more than a factor of 2 larger than at 77 = 3.2 mani- 
fested by the BRAHMS data [H . The stronger onset of 
the quantum coherence effects at 77 = 4.0 in comparison 
with 77 = 3.2 can not explain such a huge rise of nuclear 
suppression. 

Much stronger onset of nuclear effects at 77 = 4 (at 
larger Feynman xp) can be simply explained only by en- 
ergy conservation and reflects a much smaller survival 
probability of LRG in multiple quark interactions at 
larger pseudorapidity 0, [1] . 

As a demonstration of different onsets of nuclear ef- 
fects as a function of pseudorapidity we present in Fig. [2] 
predictions for nuclear suppression at different fixed val- 
ues of 77 calculating nuclear attenuation factor Rd+Auipr) 
for production of tt" mesons at y/s = 200 GeV. Changing 
the value of 77 from 3.0 to 4.0, one can see a large rise of 
nuclear suppression about a factor of 2. 

The BRAHMS Collaboration has recently reported 
new data [2j measuring nuclear effects for production of 
different species of hadrons at r] = 3.0 in d + Au colli- 
sions. These new data confirm a substantial nuclear sup- 
pression, which is similar to what has been found already 
in 2004 Using proper fragmentation functions from 
[iBl , we present in Fig. [3] the model predictions together 
with the last BRAHMS data at 77 = 3.0 for production 
of TT"*" and K'^ mesons. One can see again a reasonable 
agreement of the model with data. 
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FIG. 3; Ratio, Rd+Auipr), of identified particle production 
rates in d — Au and p — p collisions as function of pt at pseu- 
dorapidity rj — 3.0 vs. data from the BRAHMS Collaboration 
i- 



As a demonstration of the valence quark domination 
in the projectile particle leading to an enhancement in 
production of negative hadrons by deuterons, we present 
in the same Fig. [3] also the model predictions for tt" pro- 
duction at 77 = 3.0. Much smaller nuclear suppression 
clearly confirms the isospin asymmetry of leading parti- 
cle production at large forward rapidities and large px in 
d — Au collisions. 

Energy conservation applied to multiple parton rescat- 
terings in nuclear medium leads to xp scaling of nuclear 
effects [1, Q . We expect approximately the same nuclear 
effects at different energies and pseudorapidities corre- 
sponding to the same values of a;^. Such a situation is 
demonstrated in Fig. [51 where we present px dependence 
of nuclear attenuation factor Rd+Au{pT) for tt" mesons 
at different cm. energies y/s — 200, 130 and 62.4 GeV 
and such corresponding values of 77, which keep the same 
value of xp- Such a, xp scaling can be verified and inves- 
tigated in the future by the BRAHMS or STAR Collabo- 
rations in the energy and pseudorapidity range accessible 
at RHIC also for production of other species of hadrons. 



IV. SUMMARY AND CONCLUSIONS 

In the present paper we analyze a significant nuclear 
suppression in production of different species of particles 
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FIG. 4: Theoretical predictions for an approximate 
exp(77)/y^ scaling of the nuclear attenuation factor 
Rd+Au{pT) for neutral pion production in the energy and 
pseudorapidity range accessible by the BRAHMS and STAR 
Collaborations. 



at large pseudorapidities (large xp) in d — Au collisions 
investigated at present mainly by the BRAHMS P, 01 
and STAR [3] Collaborations. 

The new results of this paper are the following : 

• Using the simple formula ([6|) adopted from ref. [1| 
and based on Glauber-Gribov multiple interaction 
theory and the AGK cutting rules, we calculated 
high-pr hadron production at large x f and found a 
substantial suppression. This parameter-free calcu- 
lation agrees with recent measurements performed 
by the BRAHMS P, and STAR g Collabora- 
tions at forward rapidities in deuteron-gold colli- 
sions at RHIC. Our simple explanation is based on 
just energy conservation reflecting a small survival 
probability of LRG in multiple quark interactions 
at large xp- 

• With the same input, we explain for the first time 
very strong nuclear suppression for tt'^ production 
at ?7 = 4.0 measured recently by the STAR Collabo- 
ration 3 . This suppression is more than a factor of 
2 larger than at 77 = 3.2 (see Fig.[T|) investigated by 
the BRAHMS Collaboration [1,] for ft,^ production 



and no other models with a reasonable alternative 
description are known. 

• In order to exclude differences affected by isospin 
effects in production of different species of parti- 
cles, we performed also model predictions for nu- 
clear suppression at large pr as a function of pseu- 
dorapidity changing from 3.0 to 4.0 (see Fig. [2]). 
Predicted a huge rise of nuclear suppression with rj 
about a factor of 2 follows from much smaller sur- 
vival probability of LRG in multiple quark rescat- 
terings. 

• Using proper fragmentation functions from ref. [isj 
we performed model calculations of nuclear sup- 
pression at large pr and -q = 3.0 for 7r+ and pro- 
duction in a good agreement with the latest data 
from the BRAHMS Collaboration [2] (see Fig.[3l). 

• As a consequence of xp scaling, we predict ap- 
proximately the same nuclear effects at different 
energies and pseudorapidities corresponding to the 
same values oi xp (see Fig. 0]). Such a, xp scal- 
ing can be verified in the future by the BRAHMS 
or STAR Collaborations mainly at lower energies, 
where no effect of coherence is possible. It allows to 
exclude the saturation models or the models based 
on CGC from explanation of nuclear suppression. 

• According to xp scaling we expect similar nuclear 
effects also at midrapidities in the RHIC energy 
range. However, the corresponding values of pT for 
the produced hadrons should be much higher than 
at forward rapidities to keep the same value oi xp, 
where the nuclear suppression is studied. In this 
kinematic region, investigation of nuclear suppres- 
sion in production of different hadrons in p{d) ~ Au 
collisions is also very important because at large pt 
the data cover rather large X2 ~ 0.05 — 0.1 where 
no effect of coherence is possible. It gives another 
possibility to exclude the models based on CGC. 
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